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Thin aerofoils are prone to localized flow separation 
at their leading edge if subjected to moderate angles 
of attack a. Although Taminar separation bubbles' at 
first do not significantly alter the aerofoil performance, 
they tend to 'burst' if a is increased further or if 
perturbations acting upon the flow reach a certain 
intensity This then either leads to global flow 
separation (stall) or triggers the laminar-turbulent 
transition process within the boundary layer flow. 
This paper addresses the asymptotic analysis of 
the early stages of the latter phenomenon in the 
limit as the characteristic Reynolds number Re -> oo, 
commonly referred to as marginal separation theory. 
A new approach based on the adjoint operator method 
is presented that enables the fundamental similarity 
laws of marginal separation theory to be derived 
and the analysis to be extended to higher order. 
Special emphasis is placed on the breakdown of 
the flow description, i.e. the formation of finite-time 
singularities (a manifestation of the bursting process), 
and on its resolution being based on asymptotic 
arguments. The passage to the subsequent triple-deck 
stage is described in detail, which is a prerequisite for 
carrying out a future numerical treatment of this stage 
in a proper way. Moreover, a composite asymptotic 
model is developed in order for the inherent ill— 
posedness of the Cauchy problems associated with the 
current flow description to be resolved. 



1. Introduction 

The investigation of laminar-turbulent boundary layer 
transition is of fundamental importance in respect of the 
understanding of the complicated structure of turbulence 
and also to develop appropriate engineering models 
for the prediction of flow characteristics. Of crucial 
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theoretical as well as practical interest is the accurate calculation of lift and drag forces acting on 
aerodynamic bodies, which requires comprehensive knowledge on whether the flow is laminar 
or turbulent, attached or separated. The theory of boundary layer flows and, in particular, the 
examination of its most important issues, namely separation and transition when the Reynolds 
number is asymptotically large, have been a field of active research since Prandtl [1] presented 
his theory for laminar steady two-dimensional flows in 1904. Cornerstones of that development 
are, among others, the discovery of the singular behaviour and breakdown of the classical 
boundary layer equations near a point of vanishing skin friction (separation point in the case of 
steady flows) by Landau & Lifshitz [2] and Goldstein [3], and that of viscous-inviscid interaction 
independently made by Stewartson [4], Messiter [5] and Neiland [6] in the late 1960s, which has 
generally become known as the triple-deck theory. In conventional triple-deck problems, abrupt 
changes of boundary conditions or singular behaviour of the imposed pressure gradient initiate 
the interaction mechanism. 

On the contrary, in cases of the so-called marginal separation, a (moderate) increase in the 
smooth imposed adverse pressure gradient controlled by a characteristic parameter leads to the 
onset of the interaction process and in further consequence to localized separation. In the early 
1980s, Ruban [7,8] and Stewartson, Smith and Kaups [9] formulated a rational description of the 
local interaction mechanism now commonly referred to as the theory of marginal separation. It 
serves as the foundation of this work, which deals with the investigation of the early stages of the 
transition process triggered by the presence of laminar separation bubbles. As is well known, 
the theory of marginal separation predicts an upper bound of the control parameter for the 
existence of strictly steady, i.e. unperturbed, flows. The incorporation of unsteady effects led to 
the conclusion that the onset of the bursting process is associated either with exceeding the critical 
value of the control parameter or with the presence of a sufficient perturbation level in the case 
of below-critical conditions. Within the framework of the existing theory, vortex shedding from 
the rear of the separation bubble manifests itself in the occurrence of a finite-time singularity. 
Surprisingly, in that case recent findings strongly suggest the development of a unique blow- 
up pattern in leading order, entirely independent of the previous history of the flow [10]. The 
associated breakdown of the flow description implies the emergence of shorter scales, and the 
subsequent evolution of the flow then is described by a fully nonlinear triple-deck interaction, 
which seems to suffer finite-time breakdown as well [11,12]. The tracking of this 'breakdown 
cascade' is of particular interest and a main focus of the present investigation since it reflects the 
successive genesis of shorter spatio-temporal scales, which is a distinctive feature of the vortex 
generation process in transitional flows. 

To highlight the main issues, we restrict ourselves to the most simple case of planar 
incompressible flows. Furthermore, it is assumed that the reader is familiar with the basic concept 
of marginal separation theory, which is well established. A very detailed description can be found 
in for example [13]. Our main focus is placed on the extension of the existing theory to higher 
orders (§2) for the purpose of formulating proper initial conditions for the triple-deck stage (§3) 
and addressing the observed ill-posedness of initial value problems associated with the current 
asymptotic flow description (§4). 

2. Adjoint operator method 

In the following study, we reinvestigate the fundamental lower-deck (LD) problem of marginal 
separation theory. Here, however, we use an alternative method to that employed in the original 
studies [7-9,14] to derive the well-known similarity laws and their higher-order corrections, 
which govern the flow in the sublayer region LD (surrounded by a solid line in figure 1). To 
this end, non-dimensional quantities are introduced in the form 
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Figure 1. Asymptotic layer structure of laminar marginally separated flows including the phenomenon of bubble bursting 
(schematic). Spike formation initiated by a finite-time blow-up event at x s , modelled up to the triple-deck stage (dashed 
lines): notation and order-of-magnitude relations in terms of the perturbation parameters s = Re~ ]/2 °, o = Re~ ]/1 . Inviscid, 
irrotational upper decks (UDs), predominantly inviscid, rotational main decks (MDs), viscous boundary layer and lower deck (LD) 
regions (highlighted in grey). 



where x and y denote the coordinates in the streamwise direction and normal to the wall, t the 
time, p the pressure and xj/ the stream function. The specific flow problem under consideration is 
characterized by the dimensional reference quantities: length L, velocity Uoo, pressure poo, density 
Poo and kinematic viscosity \>oo, respectively (i.e. the unperturbed free stream values). 

Particularly, we are interested in the behaviour of the boundary layer characteristics, i.e. 
displacement thickness 8*(x) and wall shear stress r w (x), 

au =Rp -v2^ 

dy m 



5*: 
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dy m and r w : 



(2.2) 



y=0 



in the vicinity of a laminar separation bubble in the limit as the characteristic Reynolds number 
Re of the flow problem tends to infinity, 
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(2.3) 



Here, u w (x) is the velocity of the outer inviscid flow at the solid wall, y m the boundary layer 
(wall normal) coordinate and U the velocity distribution in the viscous boundary layer in the 
x-direction, which obeys the matching condition U{x,y m oo) -> u w (x). Furthermore, the usual 
scalings 
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T w 
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(2.4) 



have been used. 

The expansions of the stream function and the pressure gradient in terms of the perturbation 
parameter 

e:=Re- 1/2D ^0 (2.5) 
according to the original papers [7-9] for planar flow in the LD region in essence read 



and 
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(2-6) 



Here, the suitably scaled independent variables are denoted by 

x = s~ 4 (x — Xq), y = s~ n y and t — et, 



(2.7) 



and the origin xq of the streamwise coordinate is chosen such that it coincides with the 
point where, according to classical boundary layer theory, the wall shear stress vanishes (and 
immediately recovers downstream) as the parameter a controlling separation attains its critical 
value a c . Furthermore, pooy 3 /6 represents the separation profile, where poo > 0 is the imposed 
leading-order adverse pressure gradient; \l/ n {x,y,t) and p n {x,t) characterize the perturbation 
stream functions and the induced pressures to be determined at the levels n = l,2, ... of 
the approximation. 

Substitution of (2.6) into the Navier-Stokes equations yields to leading order 
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supplemented with the no-slip condition xj/i = dx/zi/dy = 0 at y = 0. In order to close the boundary 
value problem for (2.8), it is sufficient to require that x//i does not show exponential growth as 

y -> oo. Then yf/\ = A\y 2 /2 -\ as y -> oo follows from (2.8), where the displacement function 

A\{x,t) remains arbitrary at this stage. As a consequence, the homogeneous solution 
(eigenfunction) of (2.8) is y\r\ = A\y 2 12. The function A\ is related to the displacement thickness 
and the wall shear stress (2.2) via the expansions 



<5* ~ 55 (i 0 ) + Re~ XI \c\\A x {x, t) + q 2 x] + • • • and r w ~ Re" 1/5 ^ 3 Ai(x, t) + • 



(2.9) 



in the limit as x -> xo and a -> a c . Here, 8q(xo) represents the value of the displacement thickness 
according to classical boundary layer theory in the limit a = a c at x = xq, and the values of the 
constants q\ < 0, qi> tf3 > 0 depend on the specific flow problem under consideration. 

According to Fredholm's alternative, solutions of the inhomogeneous higher-order problems 



C^ n = b n , ft = 2, 3, 



(2.10) 



exist if and only if the right-hand sides b n are orthogonal to the eigenfunction i of the adjoint to C 
for the eigenvalue 0. This solvability condition (for n = 2) determines A\ uniquely. To be specific, 
we multiply (2.8) with the yet unknown function i(x, y) from the left and perform integration over 
the whole flow domain. Furthermore, we make use of the Fourier transform T and Parseval's 
theorem for square-integrable functions/ and g (a bar denotes the complex conjugate) 



F(f(x)) =f(k) = 
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f(x) e~ ikx dx and 



f(x)g(x) dx = 



f(k)g(k)dk. 



(2.11) 



Then multiple application of integration by parts yields 
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thus leading to the adjoint operator £ + defined by 



£+£:= 
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■ 2ikpooyi = 0, 



: 0, (2.12) 



(2.13) 



and the boundary terms 
b.t. (l,jr n ): 



y 2 d 2 f n \ didjfn d 2 l 7 



, ft = 1,2,.... 



(2.14) 



y=0 



The general solution of (2.13) may be written as 



l(k,y) = ci(k) iF 2 (l; 1/2, 3/4; Qy 4 ) + y 3 l 2 [c 2 {k)l- V ^2^Qy 2 ) + c 3 (k)I 1/4 (2^2y 2 )], (2.15) 



where pFq{a\, . . . ,a v ; b\ . . . hq,z) denotes the generalized hypergeometric function, I v (z) the 
modified Bessel function of the first kind and 



Q := 



#00 

32 



(2.16) 



From the requirement of vanishing boundary terms b.t.(£, \j/i) = 0 in (2.12), we deduce the 
homogeneous boundary condition £(k,0) = 0 from (2.14), which leads to c\{k) = 0. Furthermore, 
the suppression of exponential growth of £ as y -> oo is ensured if c$(k) = —C2(k), which actually 
results in strong decay. Using the relation K v (z) = tt[I- v (z) — I v (z)]/[2sin(v7r)] for the modified 
Bessel function of the second kind, one may write 

i = c(k)h(k,y), h := (ik) 1/8 y 3/2 K 1/4 (2Vtty 2 ). (2.17) 
Application of the procedure underlying (2.12) to higher-order problems (2.10) results in 
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Since c(k) remains undetermined, we infer the solvability condition in Fourier space to be 



d 2 h ~ dh df n 
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hb n dy = 0 



(2.18) 



(2.19) 



by using the explicit expression for the boundary term (2.14). Alternatively, inverse Fourier 
transform and application of the convolution theorem give the solvability condition in physical 
space as 
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where 
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(2.21) 



with H(x) denoting Heaviside's step function (figure 2). 

The applicability of the procedure described above, in particular in respect of the use of the 
Fourier transform, requires sufficiently strong decay of ^r n as x -> ±oo. To ensure this prerequisite, 
a shift of the form 
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(2.22) 



is performed, where a§,a\ are flow-problem-specific constants and k\ = s~ 8 (a — a c ) ~ O(l) is the 
scaled control parameter; see [7,8] for details. From the matching to the up- and downstream 

boundary layer region, it is known that A\ ^ a§\x\ +a\k\/\x\ H as ±oo; this far-field 

behaviour can be deduced from a local analysis of the classical (non-interactive) boundary layer 
equations near a point of vanishing skin friction [7]. The corresponding problem now reads 



^ 2 2 2 dx y dt 



(2.23) 



with the modified no-slip condition ^|=0, dx/z^/dy = — (A 2 — a^x 2 — 2flo0ifci)/(2poo) at y = 0. 
From (2.20) with the use of (2.21) and 
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(2.24) 




one immediately obtains — without the necessity to know the solution — * ne well-known 
fundamental equation for A\, 



A\ - ct^C - IdQdiki = - VFoo^ 
with the abbreviations 



dAl/8t -df, (2.25) 



(X-§)V4 



T(3/4) 2 3 / 4 

X = — T= and y = — -. (2.26) 

V2r(5/4) r(5/4) 

See [8,9] for the steady and [11,15] for the unsteady flow case. For the incompressible flow 
case considered here, (2.25) has to be supplemented with the upper-deck (UD) solution, i.e. the 
interaction law 



dpi Uj 



dX P0Q7T t 



d 2 A 1 /d^ 2 



df , (2.27) 



where Uoo is the velocity at the outer edge of the boundary layer evaluated at the separation point 
x = 0. 

Self-evidently, (2.19) can also be evaluated. Taking into account (2.17) and 

J y 3 / 2 + ra K 1/4 (2V% 2 ) dy = ^+ Wm y s (ikp 0 o)- {5+2m)/s r + p , m > -2, (2.28) 

one then obtains 

HA, - a 0 x - 2*0*1*1) = — T ^— j - (ifc)3/4 ? )■ (2.29) 

The Fourier transform of the pressure gradient (2.27) may be written as 

p(^\uk mA (230) 

V dX J poo 



By applying the Fourier transforms of the Weyl fractional integrals 
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(2.31) 



with the restriction 0 < v < 1, one immediately recovers (2.25). Commonly, affine transformations 
are introduced to eliminate the problem-specific constants #0/ #1/^00/ ^00 an d to underline the 
similarity law character of the fundamental equation (2.25) in combination with (2.27). With 
respect to the passage of the 'marginal separation stage' into the triple-deck stage (via finite-time 
blow-up), we prefer to keep these constants in the corresponding equations. 

Moreover, it should be emphasized that the effects of flow control elements ('smart structures'), 
such as surface-mounted obstacles and/or suction/blowing devices, can easily be incorporated 
into the analysis, for the latter, e.g. the second term in (2.20), does not vanish. The additional 
terms thus resulting may also be used for the formulation of appropriate initial value problems, 
as studied in [10] or §4. 

We are now in the position to extend the theory of marginal separation to second order. To this 
end, the explicit solution ^| of (2.23) is required, which is known in closed form in Fourier space 
only and may be found by means of a power series ansatz [13]: 
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(2.32) 



The arbitrary function b(k, t) is chosen such that ^| does not grow algebraically ~ 0(y 2 ) as 
y -> oo, and, in combination with (2.29), one can alternatively rewrite (2.32) in terms of 
generalized hypergeometric functions: 
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The far-field behaviour of and its contribution to the wall shear stress are consequently 
given by 



and 
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In order to determine the second-order correction displacement function A 2/ we follow the 
procedure described above and introduce, similar to (2.22), the shift 
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and the boundary conditions 
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Substitution of the Fourier-transformed versions of (2.36) and (2.37) into the solvability condition 
(2.19), evaluation by means of (2.32), (2.34) and (2.28), and application of the inverse Fourier 
transform yields the forced linear, fundamental equation for A 2 : 
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Here 
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are positive constants (c 2 « 0.661 009, c 3 « 0.108380). For the closure of (2.38), in addition to (2.27) 
a relationship between p 2 and A 2 is required. Investigation of the UD region (figure 1) leads to 
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see (4.15). An investigation similar to that performed in the appendix of [10] yields the far-field 
behaviour A 2 ~ OOxp 3 / 4 ) as x -> — oo and A 2 ~ 0(x 7 / 4 ) as x -> oo. 



3. Finite-time blow-up and the subsequent triple-deck stage 

As is known from [11], solutions of the fundamental equation for A\ may blow up at a finite time 
t S/ at a single point x s under certain conditions (e.g. sufficiently strong forcing for below-critical 
control parameter conditions or above-critical conditions even without any forcing, e.g. [16]). The 




specific behaviour is given by 

Ai~r- 2/3 Ai(x) + ... , pi ~r- 10/9 pi (*) + ••• and x - x s = r 4/9 x (3.1) 

as t := f s — f -> 0 + . Recent investigations show that the blow-up profile A\ or, equivalently, pi 
is unique [10] (figure 3), and consequently the following question arises. How do the initial 
conditions, effects of flow control devices, etc. (the 'history' of the flow) enter the matching 
condition to the subsequent fully nonlinear triple-deck stage first studied in [11,12] and sketched 
in figure 1? With the appropriate rescalings in the LD region 

x - x 0 - s 4 x s = o 2 X, y = (J^Y, t-s- 1 t s = aT (3.2) 

and 

4,~a 5 9(X,Y,T) + ---, j|~P00+^ + --- (3.3) 

as or := Re~ 1/7 -> 0, one obtains, after substitution into the Navier-Stokes equations, the 
fundamental LD problem 

d 2 V d 2 V d 2 V ( dV\ d 3 V 
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subject to the no-slip boundary conditions & = d^/dY = 0 at Y = 0 and the far-field behaviour 



V ~ V -^{Y + Af + 
6 



as Y -> oo and A, P -> 0 as X -> — oo. In addition, the interaction law connecting the displacement 
function A{X, T) and the induced pressure V(X, T) in the UD for incompressible flows 
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is recovered. As expected, the local displacement effect and the action of the wall shear stress (2.2) 
become intensified 
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cf. (2.9). Here c\\ < 0, q$ > 0 again denote problem-specific constants. 

Combining expansions (2.6) and the (blow-up) scalings (3.1)-(3.3) results in the 
initial/matching condition as T -> — oo, r -> 0 + : 
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Here, X= |T| 4 / 9 x, Y= \T^I 9 y, and 2i, 62 denote eigenfunctions resulting from the ansatz A2 ~ 
r _/x ^(x) for the homogeneous part of (2.38) with eigenvalues \i — (10/9, 15/9), respectively. 
Furthermore, \[f 2 represents the limiting solution ^ 2 as r 0 and, similar to (2.22), we introduce 



(3.9) 



with a yet unknown correction displacement function A 2/ which enters the asymptotic 
representation of the triple-deck displacement function 

A-p^[\T\- 6 ^ 9 A 1 + irp 10 / 9 ^ + \T\- l3 ' 9 A 2 + \T\- l5 ' 9 e 2 + • • • ] (3.10) 

as T -> —00. An asymptotic expansion of (3.4) based on (3.8) immediately leads to 
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with ^1 — 9^|/3y = 0 at y = 0. The solvability condition (2.20) then yields the equation for the 
blow-up profile 
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where pi and Ai are related via (3.6) (figure 3). Similarly, we obtain 
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Here, the indeterminate amplitudes of e\ and e 2 carry the 'history' of the flow and may be 
converted into a shift of the blow-up point x s -> x s + Ax s and £ s -> £ s + A£ s with Ax s , At s 1, 
respectively (formulation (3.4) is invariant with respect to a shift in X and T). Furthermore, the 
shapes of e\ and £2 correspond to x and £ derivatives of A\ in the limit as r 0 (figure 3) 
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To determine the blow-up profile A 2 , analogously to (2.35), ^3 is written as 
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(3.15) 



which, on further expansion of (3.4), leads to 
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13 % 

Y A 2 + -xA> 2 



-A 2 A[ 



4 9^| 4.9^1 y 9 ^2* 
+ 3~df + 9 X dxdy + 9 dy 2 ' 

supplemented with the modified no-slip boundary conditions 



y = 0: fl = ~ 



3y 



A 2 + 



3 2 Vf 2 * 



(3.16) 



(3.17) 



y=oy 



Application of the procedure that led to (2.38) here then yields 



Pi 



13A 2 + 4^ 



1/4 
PoO c i 



Pi 



(£ - £) 3 / 4 



d£+) 



-3/4 

+ Poo ^ 



jf4 



d§ 



-1/2 2A. 
oo 9 

-1/4 2C2 

.o 81 



* 3Ai + 2^ A; 



(45Ai + 50£A; + 8f 2 AiO(* " £) V4 d^ 



2c 3 



(3.18) 



Alternatively, one can derive (3.18) from (2.38) by taking the limit t t s and using the blow-up 
scalings (3.1), (3.8) and (3.10). 

With the blow-up profiles A\,Ai, the eigenfunctions e\,e2 (with arbitrary amplitudes) and 
the stream function ^| determined numerically, we are able to properly start the triple-deck 
computations. This involved task is currently under investigation. 

In general, the Cauchy problems associated with the fundamental problem (2.25) and (2.27) 
and the triple-deck stage (3.4)-(3.8) are known to be ill-posed [12,17], but can be regularized if 
the streamline curvature is taken into account; see [18-20] and the following section for a detailed 
analysis. 



4. Regularization terms 



In a recent study by the present authors [10], initial value problems based on (2.25) together with 
(2.27), i.e. 



U 2 

A\ - alx 2 - la§a\k x = — —A 



/Poo 



oo o2 



1/4 



dAx/dt 
(*-£) 1/4 



d£ +f(x,t), 



t = t 0 : A 1 (x,t)=A 10 (x), (4.1) 

were addressed numerically. Here, the additional term f(x, t) accounts for the forcing due to 
control devices such as a surface-mounted obstacle or a suction slot, respectively, and A\q{x) 
is the solution to the steady and unforced version of (4.1). Numerical solutions to the steady 
problem for different values of the parameters k\,a§,a\, poo and Uoo/ which by application of an 
affine transformation to (4.1) can be combined into a single control parameter, are presented, for 
example, in the original works [8,9]. In the analysis [10], special emphasis was placed on solutions 
that terminate in the form of finite-time singularities, and comprehensive numerical computations 
convincingly demonstrated the development of a unique blow-up structure, entirely independent 
of the particular choice of initial data. As already outlined before, further support for this 
important result was provided by an asymptotic analysis of (4.1) near the blow-up point, where 



the resulting equation admits a unique solution that is in perfect agreement with the numerical 
findings. 

On the other hand, some inconsistencies associated with the proposed Cauchy problem (4.1) 
were discovered in [15], and their repercussions on the solvability in general are discussed 
in detail in [17]. Both works mention the occurrence of instabilities in the high-wavenumber 
regime, with the latter, however, putting this in the context of an ill-posed Cauchy problem. 
More precisely, it was shown that the absolute instability of the velocity field against short- 
scale disturbances and, entailed by that, the incorrectness of the problem essentially are a direct 
consequence of the abnormal dispersion relation governing the linearized problem in the limit 
of very high wavenumbers. Using a Fourier approach in the sense of a global stability analysis 
such that 

Ai(z, t) = Aio(x) + 8 ^r-\A n (*)) + ■ ■ ■ , 8 « 1, (4.2) 
where k is real, one can derive from (2.29) and (2.26) the asymptotic result 

23/4^1/2^2 

k -> ±00 : co(k) = ^^k\k\{ik) 1 ^ + • • • . (4.3) 

This yields that, to leading order in the high-wavenumber limit, the left-hand side of (2.29) does 
not affect the dispersion relation and, above all, that the growth rates for short-scale instabilities 

k^±oo: -%(co) = D 1 \k\ 9/4 + --- , D 1 > 0, (4.4) 

are not bounded from above. Therefore, the Cauchy problem as given in (4.1) turns out to be 
ill-posed and has to be regularized. 

The essence of the above discussion is that higher-order terms in s not included so far 
must play a distinctive role in the evolution of the unsteady marginal separation process, given 
the bounded spectrum that the unsteady Navier-Stokes equations are assumed to have, but 
which the leading-order problem fails to deliver. In the light of this, the absence of short-scale 
instabilities in the numerical solutions presented in [10] has to be attributed to the applied 
(semi-)implicit scheme. Obviously, the small but, nonetheless, always present terms that result 
from the discretization and that are of the order of the truncation error proved very effective in 
regularizing the problem. For the same reason, the steady solutions serving as initial conditions 
had to be destabilized by sufficiently strong forcing in order for the blow-up to be triggered. 
On physical grounds, however, the unrealistic growth of self-excited waves should be avoided by 
analysing the Navier-Stokes equations up to higher orders rather than implementing a numerical 
scheme designed for producing the desired results. In anticipation of what follows, the solutions 
presented in [10] will nevertheless turn out to be consistent with the then well-posed problem. 

To the best of the authors' knowledge, there does not exist a rigorous asymptotic theory for the 
regularization of otherwise ill-posed Cauchy problems. The error introduced by omitting higher- 
order terms of the original Navier-Stokes equations falls beyond the scope of an asymptotic 
leading-order approach. In other words, the filtering process entailed by the temporal and spatial 
scales on which the asymptotic analysis is based may lead to the undesired result that the 
terms needed for regularizing the Cauchy problem at a certain level of approximation cannot 
be incorporated into it. They will thus enter the equations of higher order and, by that, form the 
higher-order forcing terms. Thus, the only way to regularize the original problem is to set up a 
composite asymptotic model, where terms of higher order are taken as part of the leading-order 
problem such that they can effectively contribute to stabilizing its homogeneous solution against 
short-scale disturbances. 

This method has proved to be very successful in the past (e.g. [19,21,22] and references therein) 
and will also be applied here. As revealed in these studies, for higher-order asymptotic terms 
to play the role of regularization terms, they must at least generate derivatives with respect to 
t and x of orders that are higher than those already present in the system. Inspection of the 
inhomogeneous terms in equation (2.38) for Ai when they are considered as part of the leading- 
order problem shows that, although they then do represent derivatives of higher order, essentially 



the same high-wavenumber limit (4.3), merely multiplied with a positive factor, is obtained. The 
expression for the growth rates (4.4) thus remains unchanged. However, the pressure induced 
by the second-order displacement and, to even higher order in s, by the LD stream function 
has not been examined until now, and from the arguments put forward in [21,22] it can 
be deduced that precisely the higher-order terms associated with xfa as y -> oo will give rise 
to physical phenomena that have a beneficial effect on the regularization of the leading-order 
Cauchy problem. 

In the main-deck (MD) region of the marginal separation stage, where y m = s~ 10 y, the 
expansion of the stream function resulting from the LD solutions V^i and ^2 as y -> 00, see (2.34), 
thus assumes the form 

ism{x,ym) = £ 1( Voo(ym) + ^Vml + B Y ' Vm2 + ^VmS + ^ Vm4 + £ 2 Vm5 H • (4-5) 

Substitution of (4.5) into the horizontal momentum equation then yields to leading order 



Cifmi := Voo - ^00 -fa- = -Poo + % 0 , ( 4 - 6 ) 



and additionally, 



Wm d V™ dxdy m ay^ dx 

9 2 V^ml 9^i 

>C?V /, m4 — and ^?V / "m5 — • 

9y m 3f ^ 3x 



(4.7) 



The last equation in (4.7) is interesting insofar as that here the induced pressure component 
of leading order 0(s 10 ) enters. Furthermore, in the main layer, y/x = s 6 y m /x and the pressure 
induced by V^m5 in the UD and acting in the LD forms the component of 0(s 16 ) in the pressure 
expansion. These facts suggest a pressure variation in the vertical direction for precisely this term 
of 0(s 16 ) to come into play, since its normal gradient will be comparable in size with the leading- 
order convective term in the vertical momentum equation. Consequently, the extended version of 
the horizontal pressure gradient in the MD is given by 

^ =Poo + «W + * 6 — + * 9 — + z W — +e lldj ^+ £ 12 9Pm5( *' ym '° + • • • , (4.8) 

dx dx dx dx dx dx 

where, as in (2.6), all terms representing Taylor expansions of orders higher than 0(s 4 ) have been 
omitted since they will not contribute to the results presented in the following. Most important, 
all terms in this expansion, except for /%s(x,y m , t), are identical to the corresponding components 
in the LD. As already indicated by the above order-of-magnitude argument, the asymptotic 
expansion of the vertical momentum equation then leads to the relationship 

revealing the intrusion of a normal pressure gradient. As will be shown in the following, the last 
equation in (4.7) together with (4.9) will induce a pressure response in the LD that is proportional 
to the curvature of the streamlines, i.e. the second derivative of the displacement function A\ with 
respect to the streamwise coordinate. 



After applying the rules for matching with the LD and extracting the singular parts of the 
resulting integrals, one can write the solutions to (4.6) and (4.7) as 
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and 
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dx 
dx 



dx 
dx 



00 



1 dA 1 

Poo dx 
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00^ 
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' u 2 

u 00 
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Poo 
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P00 dx ' 

■Ko 

Jo 
2 
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Al 17 " a2 ° x ) + W + PoiU ™) 



dij 



+ ^oo 

1 3A 1 
Poo dt 



^ 00 ax 



2 

PooV m 



Al 17 ~ ~ + m u ™) 



00' 



vlo^ 



U 2 

u 00 



dr] 



Vm 

IT 2 

u 00, 



(4.10) 



3 



Note that for the derivation of these expressions, the expansions y m 0 : tyoo =PooV ml & + 

(2poopoi/7\)y m H and y m -> oo : ^00 = ^ooVm H have been used (e.g. [9]). In addition, the 

equation for the normal pressure gradient (4.9) results in the expression 



ggmg = 3p5(x,t) | 1 d 3 A 1 
dx dx poo dx 3 



Woo? 



U^dri + U^yn 



(4.11) 



Here, dp$/dx denotes the induced pressure gradient of order 0(s 12 ) acting in the LD, see (2.6), 
which remains to be determined later. 

As mentioned above, the purpose of the analysis presented here is the identification of 
those pressure terms induced by the leading-order displacement function A\ that are promising 
candidates for regularizing the Cauchy problem stated in (4.1). Therefore, and for simplicity, LD 
displacement effects of orders higher in e than that represented by A 2 have been omitted in the 
MD solutions. For the same reason, furthermore, the pressure generated by the component i// m ^ 
will not be considered in the following, since neither the linear nor the quadratic nonlinear terms 
contained in ^3 will lead to derivatives of A\ with respect to x or t that are of higher order than 
that already present in the leading-order term \j/ m \ . 

The thus simplified UD expansions, with the appropriate coordinate normal to the wall now 
being y u = s~^y, are 

if U (x,yu) = £ 4 Uooy u + • 

and 



+ s u x/f ul + s 17 xlf u2 + • 



■ + s 19 x!f u4 + e 20 x!s u5 + • 



(4.12) 



U 2 

r>u(x,y u ) = —^- + ' 



1 



+ £ 10 Pu1 + £ 13 Pu2 + ■ 



+ S 15 p u4 + £ 16 pu5 + ■ 



(4.13) 



Through substitution into the momentum equations, it is easy to verify that the components 
of the stream function and the pressure satisfy 



and 



Apui = 0, 

Ap u 4 = 0, 



dy u 



+ Uoo- 



dxdt 



dx 2 



Ap u 2 = 0, 
&Pu5 = 0, 



dPu2 1T d 2 is U 2 

^7 = ^00^ 
^7 = ^00^, 



(4.14) 



where A is the Laplacian d 2 /dx 2 + d 2 /dy 2 . Relying on the matching principles and taking into 
account that apart from p m s the MD pressure terms do not depend on y m , one then obtains from 
(4.14) and (4.10) 



dpi 

dx 



POO \ dx 2 



dpi 

dx 



u2 oo n (d 2 A 2 

POO \ dx 2 



and = 

dx poo \ dxdt 



Here, H(-) denotes the Hilbert transform, i.e. 



H(f(x)) = - 



df, H 2 (f(x)) = -f(x). 



(4.15) 



(4.16) 



Consequently, the well-known leading-order result (2.27) for p\ is recovered and also the problem 
for Ai as given in (2.38) can now be stated in closed form. Furthermore, it should be noted here 
that, as a matter of course, the singular parts of the solutions (4.10) as y m —*■ oo match seamlessly 
with the corresponding counterparts in the UD. Moreover, the last equations in (4.14) and (4.10), 
respectively, together with (4.11) lead to 



-7 — (x, t, Vu = 0)- —H — V l\ = 

dx * \ dx 2 / dx poo dx 3 



U 2 00 9 3Ax 



h, 



and 



h = 



u 2 

u oo 
(^oo) 2 



4^o 2 o 
P 2 oo^ 



1 diy, h = 



(^oo) 2 
U 2 



dr], 



(4.17) 



and from (4.15) and (4.16), the LD pressure gradient dp$/dx is inferred to be given by 

dp 5 _ U 2 Q d 3 A x t 



dx poo dx 3 



ih+h). 



(4.18) 



Interestingly, results very similar to this and that for dp^/dx were found in the triple-deck analyses 
of the Blasius boundary layer with respect to the so-called lower branch stability [23] and, 
respectively, the associated acoustic radiation [24]. 

In the following, it will be assumed that t/>q 0 is smaller than Uoo throughout the MD such that 
both 1\ and I2 are positive quantities. As a consequence, the induced pressure gradient in the LD 
replacing dp\/dx from (2.27) in order to form a composite model equation when introduced into 
(2.25) can be written as 



dPc 

dx 



U: 



00 



n 



d 2 A x 



e u ^^4^C, C = (h+I 2 )>0. (4.19) 



POO \ dx 2 J poo \ dxdt J poo dx 3 

Accordingly, the modified dispersion relation derived from (2.29) for the compound model reads 

^00 7 



±00: M (k)= £ -^k\k\-^k+ , 



2 l/4„3/4 c 



leading to 



k -> ±00 : = D 2 \k\ 3/4 + 0(|/cr 1/4 ), D 2 > 0. 



(4.20) 



(4.21) 



From this, it can be concluded that the growth rates — %(co) are bounded from above and tend to 
—00 in the limit of very high wavenumbers, which, in turn, gives good reason to assume that the 
Cauchy problems based on this corrected model are well-posed. 

In order to test this assumption, we re-examined the numerical solutions presented in [10] for 
a different initial value problem by using the extended model for the pressure gradient (4.19). 
In all cases, the temporal evolution of the leading-order displacement function A\ as well as 
the development of the blow-up structure could be recalculated without difficulties. The results 



perfectly agree with those already given in that paper, even when other numerical schemes were 
applied that definitely would have failed due to the lack of inherent numerical regularization 
terms if the uncorrected model had been used. 

5. Conclusion 

This paper addresses an adjoint operator approach to the calculation of higher-order 
displacement effects in marginal separation theory In further consequence, this method is used 
to determine the initial conditions necessary for a proper formulation of the triple-deck stage, 
which is initiated due to a finite-time blow-up event in the marginal separation stage. Moreover, 
it is shown that the application of a composite asymptotic model which accounts for higher-order 
effects, such as streamline curvature, successfully leads to a regularization of the ill-posedness 
associated with initial value problems in marginal separation theory. 
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